Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 10 September 2009 (MN ]A1^X style file v2.2) 



Probing the Epoch of Reionization with Milky Way 
Satellites 



o 
o 

Oh 
C/3 



o 

U 

Oh 
6 

(N 
> 

O 

On 
O 



13 



Joseph A. Muhoz 1 *, Piero Madau 2 , Abraham Loeb 1 , and Jiirg Diemand 2 ' 3 

1 Harvard-Smithsonian Center for Astrophysics, 60 Garden St., MS 10, Cambridge, MA 02138, USA 

2 Department of Astronomy & Astrophysics, University of California Santa Cruz, 1156 High Street, Santa Cruz, CA 95064, USA 
3 Hubble Fellow 



10 September 2009 



ABSTRACT 

While the connection between high-redshift star formation and the local universe 
has recently been used to understand the observed population of faint dwarf galaxies 
in the Milky Way (MW) halo, we explore how well these nearby objects can probe 
the epoch of first light. We construct a detailed, physically motivated model for the 
MW satellites based on the state-of-the-art Via Lactea II dark-matter simulations. 
Our model incorporates molecular hydrogen (H2) cooling in low-mass systems and 
inhomogeneous photo-heating feedback during the internal reionization of our own 
galaxy. We find that the existence of MW satellites fainter than My ~ — 5 is strong 
evidence for H2 cooling in low-mass halos, while satellites with —5 > M\ > —9 were 
affected by hydrogen cooling and photoheating feedback. The age of stars in very low- 
luminosity systems and the minimum luminosity of these satellites are key predictions 
of our model. Most of the stars populating the brightest MW satellites could have 
formed after the epoch of reionization. Our models also predict a significantly larger 
dispersion in Af3oo values than observed and a number of luminous satellites with 
M300 as low as 10 6 M Q . 
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1 INTRODUCTION 

According to the standard cold dark matter (CDM) 
paradigm of cosmic structure formation, massive objects 
such as the halo of our own Milky Way (MW) grow hierar- 
chically with smaller subunits collapsing earlier and merg- 
ing to form larger and larger systems over time. Computer 
simulations have long shown that this merging is incomplete 
and that the dense cores of such progenitors may survive to- 
day as gravitationally b ound "subhalos" within their hosts 
(e.g. iMoore et al.lll999l ). Recently, state-of-the-art simula- 
tions have revealed that present-day galaxy halos are ex- 
tremely lumpy, filled with tens of thousan ds of surviving 
substructu r es on all resolved m ass scales l|Diemand et al.l 
l2008l ; ISpringel et al '"" 



2008). 



120071 . 

The predicted subhalo counts vastly exceed the number 
of known dwarf satellites of the MW, creating a "missing 
satellite problem" whose solution within the ACDM frame- 
work may lie both in the luminosity bias that affects the 
observed satellite lum inosity function ||Kqposov et al-lfeoOSl ; 
iTollerud et al.l l2008h as well as in the reduced star form- 
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ing efficiency predicted for small-mass, dwarf-sized substruc- 
ture. It is widely accepted that cosmic reionization may of- 
fer a plausible mechanism for effectively inhibiting star for- 
mation in halos that are not sufficiently massive to accrete 
warm intergalactic gas, and a number of studies have at- 
tempted to interpret the observed population of MW satel- 
lites usin g the process of early, external UV background 
feedback ([Bullock et aiB oOO: Bens on et alj2002llSomerviUel 
2002: * ' ' 

2007 
2009 



Kravtsov et al. 2004; Moore ct al. 2006: Strigari et al.l 



Simon & Gcha 2007; lMadau et al. . 2008a; B usha et all 



Maccio et ai]|2009ah . 



Lately, it has been suggested that a detailed quan- 
titative resolution of the "missing satellite problem" may 
require some "pre-reionization" suppression mechanism to 
avoid producing too many faint Galactic dwarf spheroidals 
(dSph s). Using resul t s from the Via Lactea II (VLII) simu- 
lation, iMadau et al.l (|2008bl ) showed that thousands of sur- 
viving subhalos in the halo of the MW today have pro- 
genitors massive enough for their gas to cool via excita- 
tion of H2 and fragment prior to the reionization epoch, 
which they assumed occurred around z ~ 10. In addition, 
they found that star formation in these objects must have 
been extremely inefficient converting only a small fraction of 
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their gas i nto stars or havi ng a top-heavy initial mass func- 
tion (e.g. I Abel et alj|2002 |) . Sim ilar conclusions have been 
reached bv lKoposov et all (|2009l ) . 

Inspired by these results, we develop a detailed, astro- 
physically motivated model for the formation of dSphs. We 
consider a general scenario in which H2 fragmentation can 
shut off before the suppression of atomic hydrogen cool- 
ing during reionization and include post-reionization star 
formation in the largest subhalos. This work is unique in 
that it considers the possibility that the MW was self- 
reionized from the inside out, which further suppresses the 
amount of stars produced in progenitors with M > 10 8 Mq . 
This assumption is in agreement with observations that 
the mean-free-path of ionizing radiation through interven- 
ing Lyman-limit systems (LLSs) is much less, at 2 > 7, 
than the 20Mp c distance between th e Virgo Cluster and 
the MW (|Faucher-Giguere et al.l 12009 '). Using this physical 
model with the most recent observations of the ultra-faint 
MW satellites found in the Sloan Digital Sky Survey Data 
Release 5 (SDSS DR5) as a probe of both reionization and 
early star formation physics, differentiates our work from 
previous studies that focused on the properties of dSphs. In 
this work we adopt subhalo catalogs from the one billion 
particle VLII simulation, which allow us to track the pro- 
genitors of surviving present-day MW substructure far up in 
the merger hierarchy than done before. The unprecedented 
mass resolution, combined with the fossil signatures of the 
reionization epoch in the Galactic halo, allows us to study 
gas cooling in the early universe, star formation in the first 
generation of galaxies, and the baryonic building blocks of 
today's galaxies. 

This Paper is organized as follows. In [JU W e briefly de- 
scribe the VLII simulation and develop our model. In i|3] we 
compare the results with observations, examine which ob- 
servables best constrain model components, and determine 
model parameters. Finally, we summarize our conclusions 
and discuss how they contribute to our understanding of 
reionization and high-redshift star formation in ^4] 



2 BASIC MODEL 

The recently completeted VLII simulation (Diemand et al. 
2008) uses just over a billion dissipationless particles each 
weighting 4, 100 Mq to simulate the formation of a M200 = 
1.9 x 10 12 M Q MW-sized halo in a Wilkinson Microwave 
Anisotropy Probe (WMAP) 3-year cosmology (|Spergel et al.l 
120071 ). It resolves 50,000 subhalos today within the host's 
^200 = 402 kpc (the radius enclosing an average density 
200 times the mean matter value) and tracks the merger 
history of each subhalo that survives to the present epoch 
through 400 time slices between z = 27.54 and today. In 
this study, we analyze catalogs that, for each surviving sub- 
halo, give the number of progenitors and their masses at 
a given redshift. The details o f the p rogenitor determina- 
tions are given in lMadau et alj |200 8b) . Progenitors are re- 
solved down to masses of 1O 5 M , almost a factor of 200 
be tter than the dissipa t ionles s simulation of the MW used 
bv lGnedin fc Kravtsovl (2006). This extra resolution allows 
us to incorporate physical models involving H2 cooling in 
very low-mass halos at extremely high redshift. 

Our model for assigning stars and luminosities to MW 



subhalos from the VLII catalogs involves carefully charting 
the state of cosmic hydrogen and a comparison of the cool- 
ing mass with the Jeans mass throughout cosmic time. We 
identify four epochs of star formation that can contribute to 
the stellar populations of MW dSphs: (1) stars can form at 
2 ~ 20, well before reionization, in systems large enough for 
the cooling of molecular hydrogen; (2) Hi cooling produces 
stars before reionization that are responsible for photoheat- 
ing the IGM; (3) further cooling via H I from reionization 
to 2 ~ 2 produces stars in subhalos large enough to hold 
onto their gas; and (4) stars form in the last 10 Gyr from 
metal cooling. In the following subsections, we describe the 
proces s of star formatio n in ea ch of the above epochs fol- 
lowing (Barkana~fcj]oe3 (|200ll ). and detail how we attach 
stellar populations to the simulation subhalos. In incorpo- 
rating this model, we use the set of best-fit cos mological 
param eters from the WMAP 5-year data release (iKomatsul 
120091 ). which are fully consistent with the previous results 
from WMAP3. 

2.1 Epoch of Molecular Hydrogen Cooling 

In the pristine early universe beyond z ~ 20, cooling via 
H2 is efficient in systems with total halo masses as low as 
Mh 2 ~ 5 x 10 4 Mq, while the Jeans mass is lower by an 
order of magnitude (jBarkana fc Loebl|200l| ). The resulting 
first stars create a background of Lyman- Werner photons, 
which both dissociates H2 and ac ts as positive feedback to 
replenish it l|Ricotti et alj [2002al lbh. Supernova explosions 
from these stars also begins to enrich the IGM with the 
small amount of heavy metals necessary for metal cooling. 

Ultra-faint MW dwarf satellites are a natural and un- 
avoidable consequence of H2 cooling in the pre-reionizatio n 
universe l|Ricotti fc Gnedinl 120051 ; [Bovill fc Ricottil 120081 ). 
and this mechanism can be r esponsible for the low obse rved 
metallicities of these objects l|Salvadori fc Ferrarall2008l ). We 
take advantage of the high mass resolution of VLII to trace 
the pre-reionization progenitors of these ultra-faint dwarfs, 
for the first time, down to an H2 cooling mass as small as 
10 5 Mq to make quantitative predictions about their abun- 
dance and observ able propertie s. 

Motivated bv lMadau et al.l l|2008bl ) , where an extremely 
small star formation efficiency was invoked in < 10 7 Mq 
objects at 2 = 11 to avoid overpopulating the faint end 
of the MW satellite luminosity function (LF), we allow for 
the possibil ity that H2 cooling is sup pressed earlier than 
Hi cooling l|Haiman et alj|l997t lioobT ). Earlier suppression 
could explain the very small efficiency of these small ob- 
jects in a very natural way since the halos able to cool via 
molecular hyd r ogen would be less numerous at early times. 
IWhalen et alj (2008) have suggested that this suppression 
may occur around 2 ~ 20 and result from supernova feed- 
back. 

We assume a redshift, 2h 2 , after which H2 cooling is 
quenched, and a mass threshold for halos, A/h 2 , below which 
such cooling is inefficient. We consider the VLII redshift 
snapshots at 2 = 17.9 and 2 = 23.1 as possible values for 
2h 2 and extract from the halo catalogs all surviving MW 
subhalos at 2 = with progenitors above Mu 2 a t 2h 2 . 
We further assume that a fraction /;, = Qt/^M ~ 0.16 
of the total mass in each progenitor is in gas. Due to the 
metal-poor state of the primordial gas, the first generation 
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of Population III stars has a top-heavy initial mass func- 
tion (IMF) and does not survive to the present day. Rather, 
the deaths of these stars seed the gas with traces of met- 
als and spawn a new, metal-cooled stellar population with a 
Salpeter IMF. While these metal-cooled stars are the ones 
with local relics seen today, the system initially must have 
satisfied the conditions for H2 cooling. We assume that, 
through this process, a fraction /h 2 of the initial gas is con- 
verted into stars with a Salpeter IMF; this factor accounts 
for any gas that may have been expelled by supernovae. The 
low-mass stellar population will survive to the present day 
with a visu al luminosity of My = 6 . ^ per solar mass o f 
initial stars (|Bruzual fc Charlotll2003l ; iMadau et al1l2008bh . 
For each surviving subhalo, we sum up the stars contributed 
by each progenitor at zh 2 . We assume that these stars are 
quickly incorporated into the center of the forming subhalo 
so that they are not stripped during subsequent mergers or 
tidal encounters with the MW. 



2.2 Atomic Hydrogen Cooling Before 
Reionization 

Next, we discuss the formation of stars in subhalos from H I 
cooling prior to the photoheating of the IGM at reioniza- 
tion. This cooling process is very efficient only in gas with 
temperatures above T4 = 10 4 K. Before the IGM photo- 
heats to approximately T4, this temperature requirement 
limits Hi cooling to halos whose virial temperatures ex- 
ceed T4 or, equivalently, to those with halo masses abov e 
Ma = 10 8 M Q [(1 + z)/10]" 3/2 (|Barkana fc Loebl I2001T) . 
Stars formed in this way are responsible for producing the 
A'ciump photons required on average to reionize each atom 
of Hi in the Universe. This happens when the collapse frac- 
tion of halos with M > M4 exceeds a critical threshold 
|Barkana fc Loebl 12004 ) . Thus, the condition for the aver- 
age reionization redshift of the universe, z rc i, is: 

iVclump 0. b p C rit = N-y / H I fb / CLM' M' —— (M' , Z Ici ) , 

(1) 

where dn/dM is the comoving mass function of halos, JV 7 
is the number of ionizing photons produced per baryon in 
stars, /h i is the fraction of gas turned into stars by H I cool- 
ing, and p cr it is the critical density of the universe today. 
Here, fib pcrit represents the total comoving mass density of 
baryons in the Universe to be reionized, not just those in un- 
derdense or mean-density regions. Therefore, the standard 
factor /osc, the fraction of ionizing photons that escape into 
the IGM, does n ot appear on the right- hand side of equa- 
tion CD Using the lSheth fc Tormenl (|l999h mass function for 
halos, we can solve for the combination Q = fni/N c i uulp 



1 We set a minimum stellar mass of 0.1 Mq for the Salpeter 
IMF. However, in principle, the cosmic microwave background 
creates a temperature floor to any cooling process that may in 
turn determine the m i nimum mass of stars a t high redshifts 
l|Bromm fc Loebl 1200^ : iBromm fc Larsonl 120041 1. In our discus- 
sion, we implicitly assume that stars at this lower limit would 
survive to the present day (or equivalently that this limit is be- 
low ~ IMq). In this regime, the luminosity per unit mass is not 
very sensitive to the value of the low-mass cutoff of the Salpeter 
mass function. 



given any value of z re i. For z re \ = 11, approximately the best 
fit WMAP5 value of the mean reionization redshift, we find 
Q ~ 100. 

After reionization, the IGM is rapidly photoheated to 
~ 10 4 K. Halos with masses less than the filtering mass, 
Mp, can no lo nger hold onto their ga s or accrete new bary - 
onic material l|Gnedin|[2OO0l). Whil e IKodosov et all (|2009h 
adopted the expression of iGnedinl J2000) for the remain- 
ing gas available for star formation, iBusha et al.l l)2009T j cor- 
rectly pointed out that only cold gas can actually form stars. 
The situation is made worse by infalling gas whose tem- 
perature increases by an extra order of magnitude during 
collapse. Thus, after reionization, only gaseous halos with 
virial temperatures above T5 = 10 5 K may cool via H I to 
form additional stars (see i )2.3p . 

As in £12.11 we compile a list of surviving MW subhalos 
at z — and all of their progenitors with M > M4 at z le i 
and assign to each a stellar mass M* and a luminosity (to- 
day, after more than 13 Gyr of cosmic time) of My = 6.7 
per solar mass of stars distributed with a Salpeter IMF. Ou r 
choice of IMF constrains iV 7 = 4000 ijBromm et alJl200ll ). 
As before, we assume that partial tidal stripping of their 
hosts at later times does not remove these deeply embedded 
stars. Present-day subhalos with M > Ma, progenitors at z le i 
as well as earlier progenitors at zh 2 with M > Mh 2 , as out- 
lined in $2.1\ have contributions to their stellar populations 
from both epochs. 

We now consider two alternatives for the reionization 
history of the MW galaxy -forming region (MW gfr) that af- 
fect the calculation of M* . [Alvarez et al.1 (2009) suggest that 
the MW was most likely ionized externally by radiation from 
the Virgo Cluster. However, they do not include the effect 
of intervening LLSs, which absorb external UV photons and 
may allow the MWgfr to reionize from the inside out instead. 
Here we consider and present results for both scenarios. 

If the MWgfr was reionized externally by the Virgo 
Cluster, we would expect the ionization front to cross the 
region very quickly. In this case, we can assume that the 
region reionized promptly at z Te i, and consider the subhalo 
progenitors at that redshift. The mass in stars formed by 
ea ch subhalo progenitor at 2 rG i is the same as that assumed 
bv lMadau et al 



ogemtor 
.1 (|2008bl) 



M* = /m.exAMhalo, 



(2) 



where the subscript "ex" denotes the assumption of exter- 
nal reionization. Considering only stars that formed through 
the H 1 cooling of g a s befor e 2 = 11 and assuming a Salpeter 
IMF, IMadau et all (|2008bl l found that a star formation effi- 
ciency of /h 1 ~ 0.02 for progenitors with Mhaio > 7x 10 7 Mq 
would reproduce the observed satellite LF. 

On the other hand, if the MWgfr was ionized inter- 
nally, different parts of it would have reached the critical 
collapse fraction of H i-cooling halos necessary for reioniza- 
tion at different times due to their respective overdensities 
IIBarkana fc Loebll2004l) . Thus, there should be fluctuations 
in the redshift of reionization within the MW itself about 



the mean at 



an effect that has not been considered in 



the literature. Each subhalo progenitor will photoionize and 
reheat its own gas content at the time when it has produced 
Adump ionizing photons per hydrogen atom. This condition 
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' (1 - /esc) 



f b M h 



(3) 



is met when 

M *-iv 7 (i-/ csc ) /bA4al °- 

where <?(<5) gives the ratio of N c \ uulp in the overdense pro- 
genitor to the average N c \ ump in the universe. The factor 
(1 — fesc) selects only those photons that do not escape from 
the overdense progenitor and can be used to photoionize it. 
Equation[3]tells us how much stellar mass to assign aprogen- 
itor of mass Mhaio- Each progenitor, A, has its own M > M4 
progenitors, B, that form stars at redshift z' Ici > z rc i, before 
system A collapses at z le i- Baryonic material from the re- 
gion that formed the A progenitor not converted into stars 
at z' lei is photoionized before it can be incorporated into A 
at z rc i and cannot form stars unless progenitor A eventually 
reaches a mass of M$, the mass corresponding to a virial 
temperature of T5 or a circular velocity of V5 = 52kms~ 1 . 
However, we expect most of the gas in B progenitors with 
M > A/4 at z' [ci to form stars. 

Accounting for the inside-out morphology of the MWgfr 
reionization provides a natural explan ation for the low val- 
ues of /h i. ex found by several studies (Madau et all [2008b ; 
IKoposov et al. I l2009l ; iBusha et al.l 12003 ). which assumed a 
single value for the reionization redshift of the entire MWgfr. 
Since stellar mass is a linear function of halo mass in both 
equations [5] and O both morphologies fit the same observed 
LF when 



/Hl,e 



/hi 9(5) 



fesc ) 



(4) 



If /hi<?(<5)/ (1 — /esc) is of order unity at 2 ro i = 11, we find 
that /hi, ex ~ 0.01. We can understand this physically by 
considering how, at the redshift of star formation suppres- 
sion, a lower efficiency is equivalent to a smaller amount 
of matter in collapsed objects (smaller Mhaio in Eq. [2} in 
terms of calculating how much stellar mass is produced. We 
see from this that, for a given average redshift of reionization 
in the MWgfr, the inside-out and outside-in morphologies fit 
the luminosity function equally well but result in different 
interpretations of the fitting parameters. 



2.3 Star Formation After Reionization 

We assume that cosmic reionization and the photoheating of 
the MWgfr completes rapidly after z rc i. Subsequently, only 
halos that have accumulated a mass of AI5 > 1 are able 
to hold onto their gas long enough for it to form further 
stars via Hl-cooling. We assume that this star formation 
is quenched when the subhalo begins to interact with the 
MW host. The merger histories of each surviving MW sub- 
halo at z — show only eight subhalos that achieved a 
maximum circular velocity greater than V5 at some point 
in their histories. We determine the redshift at which the 
maximum circular velocity reaches its peak. The subsequent 
reduction is due to tidal interactions during infall into the 
MW host with both nei ghboring dwarfs and the host it- 
self (|Kravtsov et al.ll2004 ) that we assume coincide with the 
quenching of star formation thro ugh ram pressure. In agree- 
ment with lKoposov et al.l l|2009l ). we find that this mass loss 
typically occurs around z = 2 — 4 (in five of our eight subha- 
los) but can happen as early as z = 8 and as late as z = 1.6. 
We add an additional mass of /s /;, M max worth of stars 



minus the mass of any stars formed in either of the first two 
epochs to each subhalo whose peak circular velocity exceeds 
V5. Here, fs is an efficiency parameter for star formation 
during this epoch that takes into account how much of the 
hot gas is cooled to form stars and how much of the gas 
mass is removed during infall. M ma x is the virial mass cor- 
responding to the maximum value of the circular velocity 
for each subhalo. We calculate the age of these new stars 
from the redshift at which the subhalo attains its maximum 
circular velocity and use this ag e to calculate the luminosit y 
of each solar mass of stars from lBruzual fc Chariot] (2003). 



2.4 Recent Star Formation 

The final episode of star formation that we include in our 
model occurs in the last 10 Gyr since z — 1.6. Since we are 
interested primarily in the physics of the early universe, we 
do not attempt to develop here a detailed model for the 
production of these stars but presu me that metal coolin g is 
involved and use the observations of lOrban et alj (|2008h to 
add additional young stars in a stochastic way. 

For each of the classical, pre-SDSS MW satellite, 
lOrban et al.l l|2008h have measured /iog , the fraction of stel- 
lar mass produced in the last 10 Gyr, and r, the mass- 
averaged stellar age. They find that the metallicities of the 
ultra-faint, SDSS satellites are so low that no star forma- 
tion is expected for these in the last 10 Gyr. As we will 
show, subhalos that achieve a maximum circular velocity 
of at least V5 are a ssociated only with the classical MW 
satellites studied by lOrban et al. (|2008l ). while stars that 
formed prior to reionization populate both the classical and 
SDSS dSphs. Therefore, we set /iog = for subhalo progen- 
itors that only formed stars before reionization to be con- 
sistent with these metallicity observations. However, setting 
this constraint only for progenitors with H2-induced cooling 
does not significantly change our results. 

The cumulative probability distribution of /iog is ap- 
proxi mately linear with a linear fit reaching unity at /iog = 
0.89 (|Orban et alJl200Sl ). Thus, for each subhalo with stars 
that formed after reionization, as outlined above, we select a 
value for /iog at random from a flat probability distribution 
within the interval [0.00, 0.89]. The mass of stars produced 
since z = 1.6 is, therefore, given by /iog/(1 — /iog) times 
the mass of stars produced from the mechanisms outlined in 
the previous three subsections. 

If we pretend that all of the stars formed in the MW 
satellites prior to z = 1.6 have an age of 14 Gyr, we can 
use the measured values of /iog and r to estimate the 
mass-weighted average age of those stars formed in the last 
10 Gyr. We find that the stars that formed after z = 1.6 
in each satellite have ages between 1 Gyr and 8 Gyr with 
an average around 5 Gyr. Since the luminosity per solar 
mass o f a stellar population doe s not vary much in this age 
range (|Bruzual fc Charlotl [iool ) . we assume a fixed age of 
5 Gyr for the recen t stars formed in each of our subhalos. 
iBruzual fc Charlotl (|2003h give the luminosity of a stellar 
population with this age as My = 5.8 per solar mass. 
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3 COMPARISON WITH OBSERVATIONS 

We are now at a position to compare the properties of the 
luminous subhalos in our model with observations of MW 
dwarf satellites. These include values for Se gue II, the newes t 
MW satellite discovered in the SDSS data [|Belokurovll2009h . 
but exclude Leo T which, at a distance of 417 kpc from 
the Sun, is outside the viral radius of the VLII simulation. 
The simulation abundance calculations represent mean val- 
ues about which there is Poisson scatter. However, for clarity 
in the plots, we instead assign Poisson errors on the obser- 
vational data as if the observed count were the mean value. 
We ignore the variations that would expected in different 
cosmological realizations of the MW and leave this analysis 
to future work. 

In H3.ll we first outline the observational biases of our 
observed sample of satellites and describe how we apply a 
similar bias to the simulated population to account for the 
completeness of the sample. We then proceed to compare the 
simulated and observed LFs of satellites in H3.2I and inves- 
tigate how this important observable constrains our model. 
Finally, we consider how well our model reproduces other ob- 
servables in £|3. 31 such as the radial distribution of satellites 
within the MW halo, the mass of a satellite within 300 pc of 
its center, and a satellite's mass-to-light ratio. 



3.1 Observational Completeness 

Our sample includes two sets of MW satellites: those dis- 
covered prior to SDSS and the ultra-faint sample found in 
the SDSS DR5. However, the observational biases for the 
two samples are not the same. Not only is the SDSS sam- 
ple complete only to a given surface brightness threshold, 
but the SDSS footprint only covers 20% of the sky. Rather 
than apply corrections to the observed distribution func- 
tions of satelli t es to account fo r the total populatio n as in 
iTollerud et ail <j2008T >. we follow lMadau et all (|2008bh in ad- 
justing our simulated population to represent those objects 
that would be detected by SDSS. 

The ob servational completen ess of the SDSS DR5 was 
modeled bv lKoposov et al l feoOSl) and lTollerud et al.l (|2008l ) 
by defining the maximum radius at which a satellite of 
a given visual magnitude would be observed. This radial 
threshold is given by: 

W = ( -4—) 10 ( - ' 6 M v-5.23)/3 MpC] (5) 

where /drs = 0.194 is the fraction of the sky covered by 
DR5. We exclude each subhalo in our simulation with a dis- 
tance beyond this threshold for its particular luminosity. 
Moreover, when plotting the luminosity and radial distribu- 
tion functions of satellites in £|3.2l and i]3.3l we correct these 
distribution by an additional factor of /drs- The observed 
distributions are also corrected so that each classical MW 
satellite contributes only /drs to the total abundance, while 
each SDSS satellite contributes fully. However, all observed 
satellites contribute equally when calculating the fractional 
Poisson errors in the distribution functions. 
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Figure 1. The luminosity function (LF) of MW satellites. Thin 
lines show abundances in the SDSS footprint with the DR5 se- 
lection threshold. The observed function is shown by the blue 
points with the non-SDSS objects each contributing a fractional 
amount /drs = 0.194 to the total. The curves show the theoreti- 
cal predictions including successively more elaborate models. The 
model given by the long-dashed, cyan curve illuminates only those 
VLII subhalos with maximum circular velocities that exceed V$ 
at some point in their histories and continue to form stars after 
reionization via atomic hydrogen cooling (with an efficiency of 
/ 6 = 0.02) and metal cooling in the last lOGyr (f 10G + 0). The 
short-dashed, red curve includes star formation in subhalo pro- 
genitors more massive than M4 at z re ; = 11.2 assuming a single 
redshift of reionization for the entire MWgfr and /m, ex = 0.02. 
The thin solid, black curve fitting the faint end of the observed 
LF additionally takes into account molecular hydrogen cooling, 
prior to suppression at zjj? = 23.1, in progenitors more massive 
than Mh 2 = 10 5 M Q with /h 2 = 0.4. The thick, solid, black 
curve represents the same model as the thin, solid one, but now 
all satellites in the MW halo out to the virial radius are show, 
not just those that meet the SDSS criteria. The long-dashed ver- 
tical line demarcates the luminosities at which pre- vs. post-SDSS 
satellites are observed. 



3.2 The Luminosity Function 

First, we would like to compare our simulated LF of MW 
satellites to the observed distribution and explore the impor- 
tance of this observable for constraining model parameters. 

The observed and simulated LFs are shown in Figure 
[U We present multiple theoretical curves, each of which in- 
cludes additional layers of the model to demonstrate the 
relationship between the model and the predictions. For ex- 
ample, taking into account only subhalos that reach a max- 
imum circular velocity of at least V5 at some point in their 
histories (see H2.3I and t|2.4|l with f$ — 0.02 roughly repro- 
duces the LF of pre-SDSS satellites. Recent, high-metallicity 
star formation in the last 10 Gyr introduces some stochas- 
ticity in the theoretical LF, but the fluctuations are smaller 
than the observed errors. However, this limited model fails 
to account for observed satellites dimmer than My ~ —9. 
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The addition of star formation in progenitors with M > M± 
before reionization at z rc i = 11.2 extends the agreement be- 
tween the theoretical and observed LFs down to My ~ —6 
when we assume /hi, ex = 0.02 in an external reionization 
model or N c \ urnp g(S) (1 — /esc) ~ 30 in an inside-out model. 
Additionally, the inclusion of pre-reionization stars lowers 
the post-reionization efficiency to f$ = 0.01. 

However, it is only with the inclusion of stars in low- 
mass molecular hydrogen cooling systems that the model 
output reproduces observations of ultra-faint satellites. This 
element of our model and its ability to explain and be con- 
strained by the data is a major new development in this work 
that distinguishes it from other studies in the literature. In 
Figure [TJ we assumed zu 2 = 23.1 and Mu 2 = 10 5 Mq with 
/h 2 — 0.4. The model not only gives the correct abundance 
of satellites brighter than the faintest known object, Segue 

I, it also provides an explanation for why no fainter objects 
have been found. To produce stars, our model requires a 
subhalo to have at least one progenitor with mass above the 
cooling mass for molecular hydrogen. The fewest stars are 
produced when the subhalo has exactly one of these progen- 
itors that is just above Mh 2 . In this case, the mass in stars 
is M* = /h 2 ft Mh 2 ~ 6600 Mq , and the luminosity today 
is My > —2.8. The existence of a satellite as dim as Segue 
1 is within the uncertainty of our approximations. While a 
precise measurement of the lower limit on the luminosity of 
MW satellites in the future will provide a tighter constraint 
on early star formation, our model implies that the exact 
value cannot be significantly below what has already been 
observed. 

We emphasize this intriguing result that star formation 
from different time periods, before the suppression of H2 
cooling and before and after reionization, are required to 
fit the theoretical LF to the data at the faint, middle, and 
bright ends, respectively. Although some objects do have 
stars from multiple epochs, they are typically dominated 
by the processes in a single part of the model. This helps 
prevent degeneracy among all of our different parameters 
and allows us to learn about star formation parameters in 
each stage almost independently of the others. 

So far, we have produced results by fixing zn 2 =23.1 
and z r ci = 11.2 and varying the efficiencies and mass thresh- 
olds to obtain agreement. However, by allowing for different 
values of the suppression redshifts, we can arrive at new sets 
of parameters that allow the model to fit the observed LF. 

Figure [5] compares the simulated LFs with reion- 
ization fixed at three consecutive VLII slices for which 
progenitor analysis is available: z rc i = 7.77, 9.14, and 

II. 2. For each of these redshifts, we find Q = 28, 46, 
and 130. We have fixed zh 2 = 23.1 and Mh 2 = 10 5 M Q 
but varied iV c iump = iV 7 fm/Q and /h 2 to find the 
values resulting in the closest fit to the data. For con- 
venience we define iV c i um p — clump g(S)/ (i- /esc) 
and / Hl = /hi g(S)/ (1 - f mc ). Our sets of fit 
parameters, then, 

{(7.77, 7, 0.05, < 0.001), (9.14, 17, 0.2, 0.1), (11.2, 30, 1, 0.4)}; 
the values of /h i, cx for an external reionization scenario can 
be calculated from equation f?] 

At earlier times, less mass in the Universe is in collapsed 
objects above the cooling mass. This results in a mono- 
tonically increasing value of Q with z le i since each baryon 
in these objects must ionize a larger number of hydrogen 
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Figure 2. The faint end of the luminosity function of MW satel- 
lites in the SDSS footprint with the DR5 selection threshold. The 
observed function is shown by the blue points with the non-SDSS 
objects each contributing a fractional amount /dR5 = 0.194 to 
the total. The curves show the theoretical predictions for different 
values of 5 rc i with zh 2 =23.1 and other model parameters op- 
timized to produce the best agreement with luminosity function 
data. The solid, black; short-dashed, red; and long-dashed, green 
curves have S rG i = 11.2, 9.14, and 7.77, respectively. The long- 
dashed vertical line demarcates the luminosities at which pre- vs. 
post-SDSS satellites are observed. 

atoms. This requires the universe to be less clumpy for fixed 
/hi- However, the reduction in the mass contained in each 
present-day subhalo's progenitors as z Te \ increases requires 
larger values of / H i to produce the same observed stellar 
mass today. This effect produces global reionization more ef- 
ficiently and allows increased clumping. Since / Hl increases 
faster with z rc i than does Q, the net result is an increasing 
JV'i 

clump 

Because later reionization allows for star formation in 
more objects near the cooling threshold before Hi cooling 
is suppressed, star formation in molecular hydrogen cooling 
halos at earlier redshifts is allowed to be less efficient to pro- 
duce roughly similar LFs. This is responsible for the reduced 
values of /h 2 that we for lower z re i. The value /h 2 < 0.001 
effectively means that the contribution of stars in molecu- 
lar hydrogen cooling systems before zn 2 to the luminosity 
of MW satellites is negligible; this results in an absence of 
satellites with luminosities fainter than My ~ —5. 

While all three considered values of z le i produced sim- 
ilar totals for the number of observed MW satellites, the 
LF slope between — 5 > My > —9 seems to favor earlier 
reionization around z rc i = 11. 

We also considered the effect of varying zh 2 and ex- 
amined two VLII slices at z = 17.9 and 23.1. We fixed the 
reionization parameters to be those best fit for z r ci = 11.2, 
but varied Mh 2 and /h 2 to find the best agreement with 
the LF data. Figure [3] shows the resulting LFs and com- 
pares them with observations. The model with parameters 
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Figure 3. Same as Figure [2] except we have fixed z rc i = 11.2 
and the solid, black and dot-dashed cyan curves show theoretical 
predictions for sets of parameters given zjj 2 = 23.1 and 17.9, 
respectively. 



set at (zn 2 , M Hj /Mq , /h 2 ) = (23.1, 10 5 , 0.4) produces al- 
most identical results as that with (17.9, 10 6 , 0.04). Allow- 
ing more time for stars to form before the suppression of Hb 
cooling is compensated for by a larger cooling mass thresh- 
old and a smaller star formation efficiency. The minimum 
cooling masses for the redshifts we considered are around 
the minimum masses typica lly assumed in simu lations of 
feedback from the first stars l|Whalen et al.ll200ct ). 



3.3 Radial Distribution, M 30 o, and M/L 

Given the degeneracies of various degrees that we have seen 
in the LF for some sets of model parameters, we attempt 
to use the other observables to learn more about the early 
universe. Figure [4] shows the radial distribution throughout 
the MW of the observed and simulated satellites for differ- 
ent pairs of fixed (z ro i,ZH 2 ) with other model parameters 
optimized to produce the best agreement with luminosity 
function data. We find that all of the zn 2 = 23.1 mod- 
els are fairly consistent with the radial data deviating non- 
negligibly only for satellites in the range 50 < r/kpc < 100 
of the MW center. The (z roi ,2 H2 ) = (11.2, 17.9) prediction 
for the radial distribution rises more sharply than expected 
for r < 100 kpc, but the statistical strength of this deviation 
is by no means overwhelming. 

We further explored whether these parameter sets were 
distinguishable in the distribution function of the maximum 
circular velocity of the subhalos. This property is difficult 
to measure observationally from the velocity dispersion of 
real satellites at a given radius from the center of the ob- 
ject. However, no significant difference was found among the 
distributions for model parameters that also fit the LF data. 

In addition to considering distribution functions, we 
would also like to test for agreement between our model 




200 300 
r/kpc 

Figure 4. The radial distribution function of MW satellites in 
the SDSS footprint with the DR5 selection threshold. The profile 
is plotted as the fraction of the total within a given radius. The 
observed profile is given by the blue points with the non-SDSS ob- 
jects each contributing an amount /drs = 0.194 to the total. The 
curves show the theoretical predictions for different pairs of fixed 
(z re j , 2jj 2 ) with other model parameters optimized to produce the 
best agreement with luminosity function data. The solid, black 
curve has (11.2, 23.1), the short-dashed, red curve (9.14, 23.1), the 
long-dashed, green curve (7.77,23.1), and the dot-dashed, cyan 
curve (11.2, 17.9). 



and observations of specific satellite properties. The mass 
enclosed within 300 pc, M300, as a function of luminosit 
has been studied us ing both obse rvational ijMateol \l99. 
iGilmore et alj [20071; Strigari et alj |2008|) and theoretica l 
|Li et al.l l200sl ; IKoposov et al.1 12OO9I ; iMaccio et al.1 l2009bh 
approaches that hinted at a characteristic mass scale for 
MW satellites of about 10 7 Mq . Both the relatively con- 
stant value of M300 and the extreme observed mass-to-light 
ratios of the faintest satellites are properties of the popula- 
tion that we would like our model to reproduce, however, we 
keep in mind that both M300 and the total mass from which 
mass-to-light ratios are calculated are difficult to determine 
observationally. 

In Figure [5] we present measu rements of M300 for MW 
satellites from lStrigari et al. I (|2008l ) and compare them with 
the output from our (z ro i, zh 2 ) = (11.2,23.1) model. Figure[6] 
similarly compares the observed mass-to-light ratios for each 
satell ite l|Mateol 1 19981 ; ISimon fc Gehal 120071 ; iMartin et ail 
l2007h with the value calculated using the simulated mass 
within its tidal radius. Both M300 and the mass-to-light ra- 
tio are plot ted as a function of s atellite V-band luminosity 
taken from iTollerud et al.l (|2008l), whose values are a little 
more up-to-date than those in IStrigari et al.l l|200Sft . While 
we continue to apply the SDSS DR5 selection criterion from 
equationO we plot all of the illuminated subhalos across the 
sky from our model in both figures rather than select only 
20 % in the SDSS field-of-view. However, we simply omit 
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known satellites for which we were unable to find reliable 
estimates of the relevant data. 

While our model predicts a fairly constant M300 as a 
function of luminosity (a variation of about an order-of- 
magnitude over six decades in luminosity) , there is still some 
disagreement at low luminosities where we find a lower av- 
erage M300 and greater scatter than observed. We find that 
values of M300 measured directly from the simulation match 
very well with those obtained assu ming, for each satellite, an 
NFW profile (|Navarro et al.ll 19971 ') fit to its simulated maxi- 
mum circular velocity and the radius at which this maximum 
is reached. A power-law fit to the resulting scatter plot of 
M300 vs. V-band luminosity, Lv, of the form M300 = /3 Ly 
yields a = 0.22 for both the simulated and NFW values of 
M3 00. This is mu c h stee per than the value of a = 0.05 for 
the lStrigari et a l. (2008) measured values of M300 with the 
up-to-date values of Lv. 

Assuming that the NFW profile is a good fit for all sub- 
halos at all times, we also calculated M300 for each subhalo 
at the redshift when the evolution in its maximum circu- 
lar velocity reaches its peak. The difference between peak 
and present values should give us an i dea of how much tida l 
stripping changes the mass estimates. iMaccio et al.l l|2009bl ) 
recently suggested that tidal stripping would result in a flat- 
tening of the M300 — Ly relation. However, our results show 
a further steepening as more low-luminosity subhalos lose 
a significant amount of their mass than do high luminosity 
ones. We find a — 0.15 for the early values of M3Q0. While it 
is true that large, low-concentration systems lose more total 
mass to tidal stripping than smaller ones do, they are actu- 
ally denser and more robust to stripping at a fixed radius. 
Additionally, if the smaller subhalos are also initially fainter, 
they must be deeper in the inner halo to be detected (i.e. 
have smaller r max ) where tidal effects are stronger, while 
larger, brighter halos can be detected even at large radii 
where tides are weak. 

Of course, an NFW profile is not always a good fit for 
all subhalos at all times. For example, if the peak in the 
maximum circular velocity is reached during a merger, one 
can get a rather large radius as well, which results in a low 
concentration and low M300. The actual mass distribution 
during a merger is quite different from NFW, and the true 
A/300 would be higher. This explains why a few of the sub- 
halos plotted in Figure [S] appear to gain mass from the peak 
until today. However, if we were able to plot the true, larger 
values of M300 at the time of the peak circular velocity for all 
subhalos, it would only show a greater steepening between 
the epoch when the peak is reached and today. 

The tension at low lumino sities between our p redic- 
tions and the measurements of IStrigari et al.l ([2008) may 
well be explained by the difficulty in determining M300 ob- 
servationally. The tracer stars in many small systems do not 
extend past 100 pc and converting the mass within this ra- 
dius to that inside 300 pc requires assumptions about the 
dark matter distribution (e.g. shape, orientation, density 
profile) and about the orbits of the stars. Additionally, in- 
terlopers and/or undifferentiated binary stars could system- 
atically skew velocity dispersion measurements, especially 
in those satellites where the dispersion is small (~ 5km/s). 
These errors could significantly inflate the mass of low-mass 
systems. 

While the observed mass-to-light ratios are calculated 
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Figure 5. Scatter plot of M300, the mass within 300 pc of the 
satellite center, versus satellite visual luminosity, Ly. Known MW 
satellites are represented by blue circles, while black x's denote il- 
luminated subhalos from our (z re i, zh 2 ) = (1 1-2, 23.1) model with 
AL300 measured directly from the simulation. Red squares show 
values of A/300 measured, for each subhalo, at the time when the 
evolution in its maximum circular velocity has reached its peak 
and assuming an NFW profile. The short-dashed blue, solid black, 
dotted green, and long-dashed lines represent power-law fits of the 
M300 — Lv relation respectively from observations, directly from 
the simulation, assuming NFW profiles for simulated subhalos 
today, and assuming NFW profiles for the subhalos when their 
maximum circular velocities peak. 



from the mass within the stellar "tidal" radius, theoretical 
values from the simulation are computed using the mass 
within the dark matter tidal radius of the subhalo. This 
means that the model predictions are upper limits on the 
observed values. This is the best we can do without modeling 
the radial distribution of the stars within subhalos. The two 
definitions of mass-to-light ratio agree only in systems where 
the stars and the dark matter are truncated by tides at 
the same radius. Despite this difference, the model roughly 
reproduces the slope and amplitude of the relation between 
satellite mass-to-light ratio and visual luminosity. We do 
appear to produce extra faint halos with ratios even larger 
than observed, which most likely is a result of the definitional 
difference between the observed and theoretical values. 

These results show that our model not only produces 
the correct abundance of luminous subhalos at each lumi- 
nosity, but it also gives subhalos with approximately the 
correct physical properties for their luminosity. That is, the 
simulated subhalos that correspond to the brightest satel- 
lites or to the ultra-faint satellites in terms of their abun- 
dance have the same physical properties as the brightest 
or faintest satellites, respectively. Our faint objects do not, 
for example, have the physical properties, such as M300 or 
mass-to-light ratio, of the brightest satellites. This consis- 
tency increases confidence in our model. 

Unfortunately, the different sets of model parameters 
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Figure 6. Scatter plot of the mass-to-light ratio of MW satellites 
versus their visual luminosities. Known MW satellites, whose ra- 
tios are calculated based on the mass within stellar "tidal" radii, 
are represented by blue circles. Black x's denote illuminated sub- 
halos from our (z re i,ZH 2 ) = (11.2,23.1) model whose mass-to- 
light ratios were calculated using the simulated mass within their 
dark matter tidal radii. 



we have been considering produce the about same range of 
M300 or mass-to- light ratio for a given luminosity. Of course, 
certain luminosities are under- or over-populated by points 
in different models, but these differences are better repre- 
sented as differences in the luminosity function as we have 
described in £ 13.21 Therefore, these scatter plots are less use- 
ful for learning about reionization and star formation in the 
early universe. 



4 DISCUSSION AND CONCLUSIONS 

In this Paper we have shown that a physical model for star 
formation applied to the subhalos of a high-resolution galac- 
tic N-body simulation can, indeed, explain the full observed 
population o f Milky Way (MW) dwarf satellites. While 
iMaccio et al.l (|2009al ) reached similar conclusions, they were 
not able to reproduce the faint end of the luminosity func- 
tion (LF) because their model did not include molecular 
hydrogen cooling in halos with masses below 10 s Mq . Our 
inside-out model of the reionization of the MW galaxy 
forming region (MWgfr) is characterized by the star for- 
mation efficiency /h 2 , cooling mass Mh 2 , and suppression 
redshift zh 2 of H2 cooling, the normalized efficiency of H I- 
cooled stars / Hl , normalized IGM clumping N^ ump , and the 
mean redshift of reionization z ro i, and the star formation 
efficiency f$ of massive systems capable of H I cooling af- 
ter reionization. Here N^ lulnp = ^ciump g{8) (1 - / csc ) and 
/h, = /hi g{S)/ (1 — /csc)- Our results are slightly at odds 
with recent measurements from the literature of M300 for 
ultra- faint systems, but we anticipate that measurement un- 
certainties may have resulted in the difference (see £|3.3p . 



In §32] we used our model to explore what the observed 
LF of MW satellites could tell us about the early universe 
and discovered that satellites of different luminosities give 
clues about star formation at different epochs: 

1. MW satellites with luminosities fainter than My ~ — 5 
can be explained by the inclusion of second-generation stars 
in low-mass halos above a cooling threshold of Mn 2 ~ 
10 5-6 Mq that were initially able to cool via molecular hy- 
drogen in the very early universe. 

We took advantage of the VLII's high resolution to probe 
this process in very small systems and allowed for the pos- 
sibility that a mechanism other than reionization was re- 
sponsible for its suppression in the early universe. We have 
shown that our theoretical subhalos illuminated in this way 
are very faint, metal-poor (due to their very early creation), 
and have extreme mass-to-light ratios and the correct abun- 
dance to be responsible for the faintest population of satel- 
lites discovered in the SDSS DR5. 

Using observations of these ultra-faint satellites to learn 
about the star formation process, we have found that molec- 
ular hydrogen cooling systems convert a fraction /h 2 of their 
initial gas into stars before star formation is suppressed 
very early in cosmic history at zn 2 . While the observed LF 
could not distinguish between models with zn 2 = 23.1 and 
zh 2 = 17.9, the radial distribution of these objects some- 
what favors earlier suppression. The best fit values of the 
cooling mass and star formation efficiency for each value 
of zh 2 , assuming cosmic reionization at z r ci = 11.2, are 
(z H2 ,M H2 /Af ,/H 2 ) = {(23.1, 10 5 , 0.4), (17.9, 10 6 , 0.04)}. 
Only if Zrci < 8 is no H2 cooling in low-mass systems re- 
quired to reproduce the observed MW satellite LF. Confi- 
dence in our models increases when we note that our values 
of the cooling threshold are consistent with what has been 
already be en assumed in simul ations of the formation of the 
first stars (|Whalen et al.|[2008h . Our results put further con- 
strains on input parameters for these studies. 

The physical difference between the two sets of 
(zn 2 , Mh 2 /Mq, /h 2 ) values that we considered may not be 
stark. The fact that fs increases for an earlier zn 2 may sug- 
gest a similarity with the inside-out reionization model of 
>i2. 21 where progenitors that reionized before the rest of the 
MWgfr have a star formation efficiency much great than 
in an in external reionization scenario. The similarities are 
intriguing, but we leave this problem to future work. 

For either value of zn 2 that we considered, ultra-faint 
satellites are clearly the sites of the very first and oldest 
in the universe, and we ag ree that searches f o r these stars 
shoul d be targeted there (|Kirbv et al ] 120081 ; iFrebel et al.l 
2009). While we do not expect stellar ages from future sur- 
veys precise enough to distinguish between formation at 
z = 20 and 2 = 8, the absence of any metal-poor stars 
older than 13 Gyr would represent a problem for our model. 

Another interesting prediction for the future is the cut- 
off that we find in the luminosity of MW satellites. We have 
shown that the molecular hydrogen cooling mass and the as- 
sociated star formation efficiency sets a minimum luminosity 
for satellites that cannot be significantly dimmer than what 
has already been observed. Interestingly, since the limit de- 
pends on the product /h 2 Mh 2 , we find the same cutoff for 
both values of zh 2 considered. While our theoretical value 
is a bit brighter than that observed to date, the difference is 
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within the uncertainties in our calculation. A precise mea- 
surement of this minimum luminosity in the future would 
be a useful test of our model, which would not be able to 
explain satellites much dimmer than Segue 1. 

2. The luminosity of satellites between — 5 > My > — 9 
is dominated by stars produced through Hi cooling of gas 
before the photoheating of the IGM at reionization. 

We considered, for the first time, a model of reionization 
in the MWgfr with an inside-out morphology where dense 
progenitors more massive than the cooling mass, Mi, reion- 
ize before the rest of the region. This prescription is based 
in the potentially large optical depth of Lyman-limit sys- 
tems between the MW and the Virgo Cluster. We showed 
that this model explains the low star formation efficiency 
found in previous studies that assumed instantaneous, ex- 
ternal reionization of the MWgfr, but both morphologies 
result in the same simulated observables. 

We used the observed LF of satellites to probe z ro i, the 
mean reionization redshift of all regions within the MWgfr. 
Depending on the choice of morphology, this represents ei- 
ther the redshift at which the MWgfr is quickly ionized by 
the Virgo Cluster or, since the overdensity of the IGM is 
small, the mean reionization redshift of the Universe. How- 
ever, if the MWgfr does self-ionize, then we can learn not 
only about the efficiency of star formation but also about 
the dumpiness of the IGM. 

After analyzing three possible values of z rc \ from 
VLII using the inside-out reionization prescription of 
i|2.2l we found corresponding sets of model parame- 
ters that best fit the LF data: (zrei, A^i ump , /hd /h 3 ) = 
{(7.77, 7, 0.05, < 0.001), (9.14, 17, 0.2, 0.1), (11.2, 30, 1, 0.4)}; 
the values of /h i, cx for an external reionization scenario can 
be calculated from equation [4] 

We find that the data is most consistent with early reion- 
ization at z rc i = 11.2, consistent with the WMAP5 central 
value, and very high values for the dumpiness of the IGM, 
star formation efficiency, and esca pe fraction. This also im- 
plies a higher star formation efficiency in molecular hydro- 
gen cooling systems at redshifts above z ~ 20 than for later 
reionization mode ls. Our result also argues for early enrich- 
ment of the IGM (|Furlanetto fc Loebll2003l '). A ssuming con- 
sisten cy with WMAP results, simulations by IWise fc Cenl 
(2009) have found that, for a "normal" (i.e. non-topheavy) 
IMF, the product /hi /esc, averaged over atomic cooling ha- 
los with virial temperatures above 8000 K, to be 0.02 with 
/esc ~ 0.4. This implies g(S) ~ 12 and that the cosmic IGM 
is relatively smooth with only a couple recombinations per 
baryon and a couple photons required, on average, to ionize 
each atom of hydrogen. 

However, it remains possible that the discrepancy between 
the data and model predictions with later reionization re- 
sults from a series of Poisson fluctuations or cosmic variance 
between VLII and the true MW history. In such a case, we 
are left with three possible points in parameter space. These 
points trace out a continuous path in the space of reioniza- 
tion parameters, but we were only able to probe discrete 
redshifts at which VLII slices have been analyzed. However, 
any external information about one of these parameters fixes 
a single point that fits the data reasonably well. For exam- 
ple, evidence that stars formed in molecular hydrogen cool- 
ing halos do not contribute to the luminosities of satellites 
today would argue for reionization around z rc i ~ 8. 



Finally, we note that the entire MWgfr need not reionize 
in one particular way; the dense progenitors could begin to 
self-ionize with radiation from the Virgo Cluster completing 
reionization for the rest of the region at some later time. The 
results we obtained by assuming an inside-out morphology 
for the entire MWgfr are valid as long as the progenitors 
of the MW satellites completely reionize themselves. This is 
because the satellite data is only sensitive to what happens 
in these objects. In this case, z r d is the mean reionization 
redshift of the Universe, since it is the time from which the 
stellar mass produced by MW progenitors to reionize them- 
selves early is calibrated. 

3. While stars formed before reionization do contribute 
somewhat to the luminosity of the brightest MW satellites 
(My < —5), this population of pre-SDSS objects can be 
explained by a combination of H I and metal-line cooling in 
the most massive subhalos after reionization. 

These satellites were represented in our model by VLII 
objects that have exceeded a maximum circular velocity of 
Vs sometime in their histories. Wh e n we included the pre- 
reionization stars and lOrban et alj (|2008l ) results to deter- 
mine the stellar mass formed in the last 10 Gyr, we deter- 
mined that a fraction fs = 0.01 of the photoheated gas re- 
tained by the largest subhalos after reionization is converted 
into stars. The majority of the gas must be kept hot or at 
low density to prevent these objects from being brighter than 
observed. 

While old stars may be found in these satellites, they are 
not prime candidates for finding the relics of the first stars. 



Although we have produced some interesting results us- 
ing MW satellite data to probe the early universe, the data 
will become sensitive to the parameters of even more compli- 
cated models if: (1) additional new satellites are discovered 
to improve the statistics of the sample; (2) the errors on 
the M300 and the mass-to-light ratio of satellites are im- 
proved; (3) the degree of cosmic variance between different 
high-resolution MW simulations is understood; and (4) new 
high-resolution simulations are developed to include baryons 
and feedback processes. Progress is already being made on 
some of these fr onts. PanSTARRS dKaiseJl2002l), the Dar k 
Energy Surv ey (iDark Energy Survey Collaboration! |2005| ). 
SkvMappe r dKellerl 120071") , a nd the Large Synoptic Survey 
Telescope l|lvezic et al.ll2008T ) will intensively probe for satel- 
lites populatin g the galactic neighborhood. Additionally, 
llshivama et al.l (|2009h have compared subhalo populations 
for many different simulations of gal actic halos and fo und 
more scatter than anticipated by ISpringel et all (|2008T ) us- 
ing a smaller sample, but their work has not yet resolved 
the progenitors that formed ultra-faint systems. Once devel- 
oped, high-resolution cosmological simulations with baryons 
and feedback will test much more specific models of reioniza- 
tion and open an avenue for comparisons with new observ- 
ables, in addition to those explored here, such as the half- 
light radius, which promises to hold interesting clues about 
the high-redshift formation physics of MW satellites. With 
these improvements, our basic methodology can be used in 
the future to further probe reionization and the process of 
star formation in the early universe. 
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